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ABSTRACT 

We present two non-parametric deprojection methods aimed at recovering the three- 
dimensional density and temperature profiles of galaxy clusters from spatially resolved ther- 
mal Sunyaev-Zeldovich (tSZ) and X-ray surface brightness maps, thus avoiding the use of 
X-ray spectroscopic data. In both methods, the cluster is assumed spherically symmetric and 
modeled with an onion-skin structure. The first method follows a direct geometrical approach, 
in which the deprojection is performed independently for the tSZ and X-ray images, and the 
resulting profiles are then combined in order to extract density and temperature. The second 
method is based on the maximization of a single joint (tSZ and X-ray) likelihood function. 
This allows to fit simultaneously the two signals by following a Monte Carlo Markov Chain 
approach. These techniques are tested against both an idealized spherical /3-model cluster 
and against a set of clusters extracted from cosmological hydrodynamical simulations with 
and without instrumental noise. In the first case, the quality of reconstruction is excellent and 
demonstrates that such methods do not suffer of any intrinsic bias. As for the application to 
simulations, we projected each cluster along the three orthogonal directions defined by the 
principal axes of the momentum of inertia tensor. This enable us to check any bias in the de- 
projection associated to the cluster elongation along the line of sight. After averaging over all 
the three projection directions, we find an overall good reconstruction, with a small 10 per 
cent) overestimate of the gas density profile. This turns into a comparable overestimate of the 
gas mass within the virial radius, which we ascribe to the presence of residual gas clumping. A 
part from this small bias the reconstruction has an intrinsic scatter of about 5 per cent, which 
is dominated by gas dumpiness. Cluster elongation along the line of sight biases the depro- 
jected temperature profile upwards at 0.2rvir and downwards at larger radii. A comparable 
bias is also found in the deprojected temperature profile. Overall, this turns into a systematic 
underestimate of the gas mass, up to 10 percent. We point out that our recovered temperature 
profiles are much closer to the mass-weighted profiles than those obtained from the X-ray 
spectroscopic-like temperature. These results confirm the potentiality of combining tSZ and 
X-ray imaging observations to the study of the thermal structure of the intracluster medium 
out to large cluster-centric distances. 

Key words: large-scale structure of Universe - galaxies; clusters; general - cosmology; mis- 
cellaneous - methods; numerical 



1 INTRODUCTION 

A precise observational characterization of the thermal structure of 
the intra-cluster medium (ICM) is of crucial relevance for at least 
two reasons. On one hand, the ICM thermodynamics is determined 
not only by the gravitational accretion of gas into the dark mat- 
ter (DM) potential wells forming clusters, but also by energy feed- 
back processes (i.e., from supernova explosions and active galactic 



nuclei), which took place during the cosmic history of the clus- 
ter assembly. On the other hand, a precise characterization of the 
temperature structure of clusters is highly relevant to infer the clus- 
ter masses, under the assumption of hydrostatic equilibrium, and, 
therefore, to calib rate clusters as precision tools for cosmol ogical 
applications (e.g., iRosati et al.ll2002l ; IVoilll2005l ; lBorganil2006. for 
reviews). 
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The study of the ICM properties has been tackled so far 
through X-ray observations. Data from the Chandra and XMM- 
Newton satellites are providing precise measures of the temperature 
and surface brightness profiles for a fairly large number of nearby 
(2^ 0.3) clusters, reaching z ~ 0.5 for the brightest objects (e.g. , 
Piffaretti e t al."2005': 'Pra tt & Amauj|2005l : IVikhlinin et al.ll2005l : 
Kotov & V ikhlinin 2006). 

These observations have indeed allowed to trace in detail the 
mass distribution in galaxy clusters for the first time. However in 
the X-rays the accessible dynamic range is limited by the p^^j^ de- 
pendence of the emissivity which causes measurements of the tem- 
perature profiles to be generally limited to 2-3 core radii, extend- 
ing out to rsoo only in the most favorable cases. This is not the 
case f or clusters' studies performe d with the Sunyae v-Zeldovich 
effect ( Sunvaev & Zeldovich '1972", tSZ hereafter; see lSirkinshawl 
[l999; Carlstrom et al. 2002 for reviews). Since the tSZ signal has a 
weaker dependence on the local gas density, it is in principle better 
suited to sample the outer cluster's regions, which can be accessed 
by X-ray telescopes only with long exposures and a careful charac- 
terization of the background noise. Clusters are currently observed 
through their thermal SZ (tSZ) signal and tSZ surveys of fairly large 
area of the sky promise to discover in the next future a large number 
of distant clusters out to z'SL 1. 

Thanks to the different dependence of the tSZ and X-ray emis- 
sion on the electron number density n^, and temperature Tg, the 
combination of these two observations offers in principle an alter- 
native route to X-ray spectroscopy for the study of the structural 
properties of the ICM. Indeed, while the X-ray emissivity scales 
as nj.K{T) (where A(r) is the cooling function), the tSZ signal is 
proportional to the gas pressure, neTe, integrated along the line- 
of-sight. Recovering the temperature structure of galaxy clusters 
through the combination of X-ray and tSZ data has several ad- 
vantages with respect to the more traditional X-ray spectroscopy. 
First of all, surface brightness profiles can be recovered with a lim- 
ited number (~ 10'^) of photons, while temperature profiles require 
at least ten times more counts. Therefore, the combination of X- 
ray surface brightness and tSZ data should allow to probe more 
easily the regimes of low X-ray surface brightness (i.e. external 
cluster regions and high-redshift galaxy clusters), which are hardly 
accessible to spatially resolved X-ray spectroscopy. Furthermore, 
fitting X-ray spectra with a single temperature model is known to 
provide a temperature estimate which is generally biased low by 
the presence of relatively cold clumps embedded i n the hot ICM 
atmosphere jMazzotta et ai]|2004 IVikhlinin|[200^ . On the other 
hand, combining X-ray and tSZ does not require any spectral fit- 
ting procedure and, therefore, yields a temperature which is basi- 
cally mass-weighted. 

The combination of X-ray and tSZ observations is cur- 
rently used_to^ estimate the angular diameter distance of clusters 
(e.g., iBoname nte et al. 2006; Ameglio et al. 200^. , and references 
therei n) and to recover the gas mass fraction (e.g.. lLaRoaue et al] 
I2OO6I) . Clearly, performing a spatially-resolved reconstruction the 
thermal structure of the ICM requires the availability of high- 
resolution tSZ observations with a sub-arcmin beam size, with 
a sensitivity of few /iK on the beam. Although observations of 
this type can not be easily carried out with millimetric and sub- 
millimetric telescopes of the present generation, they are certainly 
within the reach of forthcoming and planned instruments of the next 
generation, based both on interferometric arrays (ALMA: Atacama 



Large Millimeter Arra}^) and on single dishes with large bolometer 
arrays (Comell-Caltech Atacama Telescope: CCATo; Large Mil- 
limeter Array: LMlQ). 

Combining X-ray and tSZ data to reconstruct the three di- 
mensional gas density and temperature structure of galaxy clusters 
is not a n ew idea and di f ferent authors have proposed different ap- 
proaches. IZaroubi et al.l ( Il998l) used a deprojection method, based 
on Fourier transforming tSZ, X-ray and lensing images, under the 
assumption of axial symmetry of the cluster. After applying this 
method to simple analytical cluster models, they concluded that 
the combination of the three maps allows one to measure indepen- 
dently the Hubble constant Hq and the inclination angle. This same 
method was then applied bv lZaroubierZI ( l200lf) to cosmological 
hydrodynamical simulations of galaxy clusters, who found that a 
reliable determination o f the cluster bary on fraction, independent 
of the inclination angle. iReblinskvl J200ol) applied a method based 
on the Richardson-Lucy deconvolution to combined tS Z, X-ray 
and w eak lensing data to a set of simulated clusters. lDor6 et al.l 
(I2OOI') used a perturbative approach to describe the three dimen- 
sional structure of the cluster, to combine tSZ and lensing images. 
In this way, they were able to predict the resulting X-ray surface 
brightness. After testing their method against numerical simula- 
tions of clusters, they conclude that t he DM and gas d istributions 
can both be recovered quite precisely. iLee & Sutol ^2004) proposed 
a method, based on assuming a polytropic equation of state for gas 
in hydrostatic equilibrium, which allowed them to recover the three 
dimensional pro files of clusters usin g the tSZ and the X-ray sig- 
nals. jP uchwein & BartelmannI j2006l) applied the same method of 
iReblinskvi t200Ci) to deproject X-ray and tSZ maps, so as to recover 
the gas density and the tempe rature structure of cluste rs, under the 
assumption of axial symmetry. Cavali ere & Lapil fe006) applied the 
combination of tSZ and X-ray observations to recover determine 
the ICM entropy profile. 

As for applications to real clusters, [zhang & Wul j2000l ) com- 
bined X-ray surface brightness and tSZ data, for a compilation of 
clusters, to estimate the central cluster temperature, and found it to 
be in r easonable agreement with t he X-ray spectroscopic determi- 
nation. |Pointecoutem|et^] |200^ used ROSAT-HRI imaging data 
of a relatively distant cluster {z ~ 0.42) with tSZ observations to 
infer the global temperature of the system. 

De F ilippis et alj j200A combined X-ray and tSZ data to 
constrain the intrinsic shapes of a set of 25 clusters. By ap- 
plying a deprojection method bas ed on assuming the /3-model 
dCavaliere & Fusco-Femiaiiol Il976h . they confirmed a marginal 
preference for the clusters to be aligned along the line-of-sight, thus 
conclu ding tha t X-ra y selection may be affected by an orientation 
effect. H erend ( l2007h analyzed the potentiality of combining tSZ, 
X-ray and lensing data to constrain the 3D structure of the clusters. 
He found that these data are enough to determine the elongation 
along the line of sight (together with the distance), without how- 
ever fully constraining shape and orientation. 

Some of the detailed methods applied to numerical cluster 
models account for the presence of a realistic noise in the tSZ and 
X-ray maps. However, they generally do not present any detailed 
assessment of how this noise determines the uncertainties in the de- 
projected profiles, which ultimately characterize the ICM thermo- 
dynamics. Having a good control on the errors is especially crucial 



^ 'http : / / www . eso . org/pro ject s / alma/ 

^ http://astrosun2.astro.corneIl.edu/research/projects/atacama/ 
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in any deprojection technique, since errors at a given projected sep- 
aration affect the deprojected signal in the inner regions, thereby 
introducing a non-negligible covariance in the reconstruction of 
the three-dimensional profiles. 

In this paper we discuss a method to recover the three- 
dimensional temperature and gas density profiles from the joint 
deprojection of X~ray surface brightness and spatially resolved tSZ 
data, testing its performance against idealized spherical clusters 
and full cosmological hydrodynamical simulations. This method is 
based on the assumption of spherical symmetry, but do not assume 
any specific model for the gas density and temperature profiles. 
We will describe two different implementations. The first one is 
analogous to that alre ady applied to deproject spectroscopic X-ray 
data (e.g., iKriss et aU1983) and is based on assuming a onion-like 
structure of the cluster, in which projected data of X-ray and 
tSZ "fluxes" are used to recover gas density and temperature in 
the external layers and then propagated to the internal layers in a 
iterative way. The second implementation is based instead on a 
multi-parametric fitting procedure, in which the fitting parameters 
are the values of gas density and temperature within different 
three-dimensional radial bins. The values of these parameters are 
then obtained through a Monte Carlo Markov Chain maximum 
likelihood fitting by comparing the resulting projected X-ray and 
tSZ profiles to those obtained from the maps. As we shall discuss in 
detail, this second method naturally provides the error correlation 
matrix, which fully accounts for the covariance between error 
estimates at different radii and among different (i.e. gas density and 
temperature) profiles. The quality of the X-ray data required by 
our methods are basically already available with the current gener- 
ation of X-ray telescopes. As for the tSZ data, exploiting the full 
potentiality of the deprojection requires spatially resolved data. For 
illustrative purposes, we will assum e the forecast observing con- 
ditions and sensitivity of the CCAT jSebring et alj|200Q . see also 
rhttp://wwwsubmm.caltech.edu/~sradford/ccat/doc/2006-01-ccat:feSiibilB^i'l_?l^ 
although our computations can be easily repeated for other tele- 
scopes. 

The plan of the paper is as follows. In Section 2 we describe 
the two implementations of the deprojection method, while we de- 
scribe in Section 3 their application on a spherical polytropic f3- 
model. Section 4 presents the results of the analysis on the hy- 
drodynamical simulations of clusters. The main conclusions of our 
analysis are summarized in Section 5. 
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Figure 1. Illustration of the onion-skin model adopted for the geo meti'ical 
deprojection (see text in Sect. 2.1; adapted from lMcLaughlinll 19991) . 



the deprojection method within the i-th shell. In our analysis will 
be proportional to either UeTe for the tSZ signal, or to n^A{T) for 
the X-ray emissivity. In this way, the contribution of the i-th shell 
to the surface brightnes jf| in the ring j of the image will be given 
by Si^j = e; • Vi.jjAj, where the matrix Vij has as entries the 



values of the volume of the shell i which is projected on the ring j, 
whose area is Aj. By definition, s^j = for j > i. Accordingly, 
the surface brightness S'j in the ring j can be obtained by summing 
up the contributions from all the shells. 



2 THE METHODS OF DEPROJECTION 

2.1 The geometrical deprojection technique 

The first method that we apply to recover the three-dimensional 
profiles of temperature and gas density is based o n a ge omet- 
rical technique originally introd uced bv iKriss et al.l jl983), and 
subsequently adopted by (e.g., ISuotd [2OO0I : Ettori et alj 20021 : 
[Morandi et al. 2007) to deproject X-ray images and spectra of 
galaxy clusters. This method of geometrical deprojection is fully 
non-parametric and allows to reconstruct the 3-dimensional pro- 
file of a given quantity from its 2-dimensional observed projection, 
under the assumption of spherical symmetry. 

Following Kriss et al. (1983), the cluster is assumed to have 
a onion-like structure (see Figure [T](, with A'^ concentric spherical 
shells, each characterized by uniform gas density and temperature 
within it. Therefore, the cluster image in projection is divided into 
rings, which are generally assumed to have the same radii of the 3D 
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where the sums extend over the radial bins. The deprojection 
amounts to invert the above equation, i.e. to recover the values of 
from the observed projected signal S'j. We refer to Figure[T]to il- 
lustrate how this deprojection is performed in practice. Let the shell 
i, limited by r-i and r^+i, be the outermost one. Then, from the sur- 
face brightness S'i in the ring i (limited by Ri and -R,;+i), one can 
directly compute the emissivity of the shell i simply by knowing 
the volume of the region (a) and the area of the ring. In this case, 
the sum in eq.[T]has only the term j = i — N . The adjacent inner 
ring, having index i — 1 and limited by and Ri, takes instead 
a contribution from both the i — 1 and i shells. The former is com- 



^ For the sake of clarity, we indicate here with surface brightness the pro- 
jected quantity, which can be both a genuine X-ray surface brightness and 



the tSZ signal. 
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puted by multiplying the emissivity of that shell by the volume of 
the region (b). After subtracting it, the only remaining contribution 
is that of the sell z — 1 from which the emissivity e^-i is computed. 
This procedure is then repeated from ring to ring down to the center 
of the cluster. 

For this simple scheme to be applied, one requires to have im- 
ages extended out to the true external edge of the cluster, i.e. out 
to the radius where the surface brightness goes virtually to zero. 
Clearly, this situation is never attained in practical applications for 
at least two reasons. First, clusters are always embedded in a large- 
scale cosmic web, which makes it difficult to define a sharp outer 
boundary. Second, and more important, both instrumental and cos- 
mic backgrounds often dominate the genuine signal from the clus- 
ter well before its virial boundary is reached. 

To overcome this problem, it is then necessary to take into ac- 
count the emission from the gas, which extends outside the A'^-th 
shell. This emission does not have a corresponding ring in the im- 
age but can give a non-negligible contribution to the surface bright- 
ness in all r ings. To account fo r this contribution, we follow the 
approach of iMcLaughlinI ( 1 1999 1, who modeled the volume emis- 
sion from the gas beyond the last observable annulus as a power 
law, e(r) oc r^" (we refer to their Appendix A for a more detailed 
description). The idea behind this method is that the exact contri- 
bution to each ring from the external part can be calculated by in- 
tegrating the volume emission e(r). Then, the normalization of the 
power law shape of e(r) is fixed by the requirement of matching 
the total surface brightness of the last ring. This correction can be 
expressed as an additional term to eq.(|TJ, which is proportional to 
the surface brightness of the last ring: 

Sj =S',+ f, ■ Sn , (2) 

Here, fj is a geometrical factor which is uniquely specified by the 
values of the limiting radii of the j-th ring and by the exponent a. 
Eq.^ must be actually interpreted as a set of 27V equations, which 
corresponds to the separate deprojection of the tSZ and of the X- 
ray signal, each performed for A'^ radial bins. The geometrical de- 
projection is then performed by inverting each set of A'^ equations, 
starting from the outermost bin and proceeding inward. This pro- 
cedure provides the radial profiles of UeTe and of n'^A{T), whose 
combination finally gives the 3D profiles of electron number den- 
sity and of temperature. We emphasize that the temperature so ob- 
tained is the actual electron temperature and not the deprojected 
spectroscopic temperature, usually obtained from the fitting of X- 
ray spectra. 

Given the iterative nature of this procedure, the uncertainty 
associated to each ring propagates not only to the corresponding 
3D shell, but also to all the inner shells. For this reason, it is very 
difficult with this method to have a rigorous derivation of the sta- 
tistical uncertainties associated to the deprojected profiles. This is 
particularly true for the X-ray profiles, that also involve a deriva- 
tive of the cooling function with respect to the temperature. The 
commonly adopted solution is based on realizi ng MonteCarlo sim - 
ulations, over which to compute the errors (•e.g. lEttorietalj2002l) . 

Furthermore, errors associated to different radial bins are not 
independent. This is due to the fact that the projected signal in a 
given ring is contributed by several shells. The resulting covariance 
in the signals recovered in different shells is not provided by this 
deprojection method. This is a rather important point on which we 
will come back in Section 4. 



2.2 The maximum likelihood deprojection 

This technique is based on performing the deprojection by maxi- 
mizing a likelihood function, which is computed by comparing the 
observed tSZ and X-ray profiles with the ones obtained by project- 
ing the onion-skin model in the plane of the sky. This approach 
offers more than one advantage with respect to the geometrical de- 
projection, described in the previous section. First, the deprojec- 
tion of both X-ray and tSZ profiles is performed simultaneously, 
directly obtaining the whole density and temperature profiles and 
their errors. Second, besides the variance, it is also possible to com- 
pute the correlation matrix for all parameters, without any extra 
computational cost. Finally, it is possible to introduce in the likeli- 
hood extra terms in order to improve the accuracy and robustness 
of the technique. As we shall describe in the following, we adopt 
a regularization technique, which is based on imposing a suitable 
constraint to the likelihood function, to smooth out spurious oscil- 
lations in the recovered profiles induced by the covariance in the 
parameter estimate. 

The definition of the likelihood is the most important part of 
the whole procedure. We define a joint likelihood for the tSZ pro- 
file, Ctsz, and for the X-ray surface brightness profile, Cxray, 
also including a term associated to the regularization constraint, 
Creg- Since these three terms are independent, the total likelihood 
is given by the product of the individual ones: 

C = Ctsz ■ Cxray ' .C^eg- (3) 

For both the tSZ and the X-ray profiles, the take the Gaussian form 
for the likelihood, 

ln{CtSZ.X-ray) ^ =~2^\ 7 ) ' 

i 

where Oi are the values of the profiles, in the i-th bin, measured 
from the maps, while Mi{x) are the model-predicted profile val- 
ues, as obtained for the set x of parameters. Finally, at is the un- 
certainty on the measured values Oi . 

While the Gaussian expression is adequate for the tSZ signal, 
its application for the X-ray, instead of the Poisson distribution, re- 
quires the number of photons sampling the surface brightness map 
in each radial shell to be large enough to neglect the Poisson noise. 
As we shall discuss in the following, even in the outermost rings, 
we always have at least 20 photons in the "noisy" X-ray maps. 

For the regularization con straint, we ado pt the Philips- 
Towmey regularization method teouchetl [1995I and references 
therein ). This method has been already used also bv lCroston et al.l 
( l2006h to deproject X-ray imaging and spectral data. The method 
consists of minimizing the sum of the squares of the fc"* -order 
derivatives around each data-point, so as to smooth out oscilla- 
tions in the profiles. Here we choose to minimize the second-order 
derivative, since we aim to eliminate fluctuations in the profiles, 
but not the overall gradient. As we shall discuss in the following, 
such oscillations are due either to genuine substructures or to noise 
which propagates from adjacent bins in the deprojection. The local 
derivative of the function Xi at the i-th radial interval is computed 
by fitting it value and the values at the adjacent points, Xi-i and 
Xi+i), with a second order polynomial. Let Vi be the value of the 
equally-spaced cluster-centric distances, at which the profiles are 
sampled, and Ar the spacing. Then, the regularization likelihood 
can be cast in the form 

M^.eJ = ^ ( ^1 j = 

i=3 
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The quantity between parenthesis in the line formula is the exact 
value of the local second-order derivative around r^. All the con- 
stant factors are included in the coefficient A, which is called the 
regularity parameter. The choice of its value is determined by the 
compromise one wants to achieve between the fidelity to the data 
(low A) and the regularity of the solution (high A). A small A value 
will give an inefficient regularization, while a too high A will force 
the profile to a straight line, especially if the signal-to-noise ratio, 
S/N, is low. We apply the regularization constraint only to the tem- 
perature profile, which is that generally showing large oscillations, 
while the density profile has always a rather smooth shape. The sum 
in eq.® starts from i = ?> since we prefer to exclude the innermost 
point from the regularization procedure. 

With this approach, the values of the 3D gas density and tem- 
perature profiles are computed at A'^ = 15 radii each. Therefore, 
the total number of parameters to be determined with the maximum 
likelihood approach is 30. In order to optimize the sampling of such 
a large parameter space, we adopt a Ma rkov Chain Monte Carlo 
(MCMC) fitting technique tNeal 1993. : Gilks et alJI 199^ : iMacKavl 
[1993) . 

What the MCMC computes is the (marginalized) distribution 
of each parameter of a set, Xi (i.e. the values of density and tem- 
perature into each bin), for which the global (posterior) probability 
P{x), which is proportional to the likelihood function, is known 
at any point in the parameter space. In the case of an high number 
of parameters or of a particularly complex P{x), this is quite diffi- 
cult to be done analytically, or simply computationally very expen- 
sive. Instead, the MCMC performs the exploration of the parame- 
ter space with a limited computational cost, thanks to an iterative 
Monte Carlo approach, by sampling the P{x) distribution. At each 
iteration, new values of the parameters are drawn from a symmetric 
proposal distribution, that in our case is a Gaussian, 



q[Xi,Xi) oc e ^ ' *' ' ' 



(6) 



Here Xi and Xi are the entries of two vectors, having 30 components 
each, which represent the updated and the old values of the fitting 
parameters, respectively. The parameter at determines the possible 
range for Xi given x . 

After the likelihood function is computed for a new set of 
parameters x, these new values are accepted or rejected with 
a probability (A) give n by the so-called Metropolis criterion 
( [Metropolis etal.lll953h : 



A{x, x) 



P{x) 



(7) 



where P{x) is the distribution sampled by the MCMc[fl 

The width of the proposal distribution appearing in eq.(|6ll, a^, 
determines the behavior of the chain: a small value of ai increases 
the acceptance rate since the new proposed value is close to the 



5 iHastingj h97d) has generalized this treatment to non-symmetric pro- 
posal distributions, by adding a factor in equation [7] which takes into ac- 
count the proposal distribution q{x, x): 



A(x,x) = rain { 1, — ' ' } , 
^ ' I ' P{x) q{x,x) \ ' 

This is called the Metropolis-Hastings criterion. 
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old one, while a high value provides a faster exploration of the pa- 
rameter space. Our criterion to choose the values of ai is that the 
resulting acceptance rate, given by eq.Q, is around 10 per cent. 

Each parameter is allowed to vary within a finite interval, in 
order to avoid that the MCMC finds secondary maxima in unphysi- 
cal regions of the parameter space. As for gas density, we allow it to 
vary within a large range, 0.1 < rie < 10~® cm~"^. Since density 
is mostly constrained by the X-ray signal, which is proportional to 
iT^, it is always fairly well constrained and the above large interval 
of variation does not create convergence problems in any of our ob- 
jects. The upper and lower limits allowed for the temperature are 
25 keV (never reached along the chain) and 0.5 keV. Even though 
none of our clusters reach such low temperatures within the virial 
radius, the exploration of the parameter space during the Markov 
Chain run could reach such a low temperature regime. When this 
happens, the rapid drop of the cooling function A(T) below 0.5 
keV generates a maximum in the likelihood probability distribu- 
tion, with unphysically low temperature and very high density. 

The iterative procedure described above is repeated until a 
suitable number of new sets of parameters are accepted in the chain 
(typically ^ 5 x lO**). In this condition, the frequency of the oc- 
currence in the chain of the i-th parameter Xi approaches its true 
probability distribution, P{xi). Note that each parameter distribu- 
tion is already marginalized over the distributions of all the other 
parameters. 

We perform the statistical analysis of the chain by using the 
code getdist of the COSMOMC package l lLewis & Bridlel2002l) . In 
addition to a complete statistical analysis of the chain, the code per- 
form s a series of convergenc e tests: the Gelman & Rubin R statis- 
tics I Gelman & Rubir]|l992h . the Raftery & Lewis test teaftervl 
'2003') and a split-test (which essentially consists in splitting the 
chain into 2, 3 or 4 parts and comparing the difference in the pa- 
rameter quantiles). We check the convergence of our result against 
all these three tests. 



3 APPLICATION TO AN IDEALIZED CLUSTER MODEL 

In order to investigate the presence of possible systematics in the 
geometrical deprojection technique, we carry out a test on an ideal 
cluster model. We construct this model cluster by assuming the (3~ 
model for the gas density profile, with an effective polytropic equa- 
tion of state to define the temperature profile: 



To 



-3/3/2 



VneO 



7-1 



-3/3(7-l)/2 



(9) 



The values of the model parameters are fixed as follows: f3 = 0.8, 
7 = 1.2 for the effective polytropic index, n^o ~ 3 ■ 10^^ cm~^ 
for central electron number density. To = 4 keV for the central 
temperature, Vc = 200 kpc for the core radius. The "virial" ra- 
dius, which represents here the largest cluster-centric distance out 
to which the profiles are followed, is fixed at r^ir ~ 2 Mpc. We 
assume the cluster to be placed at redshift z = 0.1. 

We create tSZ and X-ray images of this model in X-ray and 
tSZ. The maps are composed by a grid 512x512 pixels. The phys- 
ical dimension of the pixel is ~ 16 kpc, which corresponds to 
an angular scale of ~ 4 arcsec at this redshift. In realizing the 
map of X-ray surface brightness, we adopt for the cooling function 
the pure bremsstrahlung expression A(r) = Ao(T/ro)°-^ where 
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Figure 2. Results of the geometrical deprojection of the analytical model: density profile (left panel), temperature profile (right panel). The dashed line 
represents the true profile. The circles (triangles) connected by a solid line represent the reconstructed profiles without (with) the inclusion of the noise. 
Errorbars represent la deviations over 1000 Monte Carlo resamplings over the data (see text for details, Section lTTl 



Ao = 5 ■ 10^"^^ erg/s cm'^. We prefer this simple formula since at 
this stage our main interest is to investigate the systematics of the 
deprojection itself, rather than the uncertainties introduced by the 
dependence of the X-ray emissivity on the temperature (which we 
expect anyway to be quite small). 

We adopt the same binning strategy for both the ideal clus- 
ter and for the simulated objects, that we shall describe in the next 
Section. The first bin is taken from r = to r = 0.05r„ir which 
always corresponds to ^ 100 kpc in our set of simulated clusters. 
Then, we compute the profile in 10 (15) bins out to 7?5oo (Rvir) 
which are equally spaced in logarithm. This choice represents a 
good compromise between the needs of accurately resolving the 
profile and of having an adequate signal-to-noise (S/N > 5) in each 
bin. We point out that a proper binning criterion is important in or- 
der to get an unbiased reconstruction of the profiles. One should 
keep in mind that the spherical shells are assumed to have homoge- 
neous gas density and temperature structures, thus neglecting any 
internal radial gradient. On the other hand, the portion of each shell, 
which is projected on the corresponding ring in the image, is lo- 
cated at a larger radial distance from the center than the portion of 
the same shell which is projected into the inner rings. Therefore, if 
the bin width is comparable or larger larger than the scale length of 
the internal radial gradient of the shell, the emissivity contributed to 
the correspondent ring is lower than expected from a homogeneous 
gas density, while it is larger for all the inner rings. As a conse- 
quence, the emissivity of the shell is underestimated while that of 
all the inner shells is overestimated in order to correctly reproduce 
the cluster image. 

In order to check for the presence of such systematics, we first 
apply the deprojection technique in the case of an ideal observa- 
tion, free of any noise. The reconstructed density and temperature 
profiles are shown in Figure|2] The reconstruction in this extremely 
idealized case is excellent, with very small or no deviations in all 
bins. Larger deviations are in the outermost bins and are related 
to the subtraction of the contribution of the fore-background con- 
taminations. This contamination is due to the fact that the /3-model 



used to produce our maps ideally extends out to infinity. Neverthe- 
less, all deviations are smaller than a few percent and are negligible 
with respect to any observational noise. The results obtained in this 
test case show that taking equally log-spaced bins is in fact a good 
choice. 



3.1 Geometrical deprojection of the noisy maps 

The case of noiseless observations discussed in the previous sec- 
tion is highly idealized. The impact of including a realistic noise is 
instead very important and cannot be neglected. The recipe to add 
noise to the maps, that we describe here, will also be used in the 
study of the simulated clusters, discussed in Section 4. 

As already mentioned in the Introduction, recovering detailed 
temperature profiles from the combination of X-ray imaging and 
tSZ data requires both of them to have an adequate spatial reso- 
lution. While this is certainly the case for the present generation 
of X-ray satellites, the combination of good sensitivity and spa- 
tial resolution for tSZ observations should await the next genera- 
tion of sub-millimetric telescopes. For the purpose of our analysis, 
we model the noise in the tSZ maps by using as a reference the 
performances expected for the planned Cornell-Caltech Atacama 
Telescope (CCAT), which is expected to start operating at the be- 
ginning of the next decade lSebring et alj2006l) . The telescope will 
be a single-dish with 25 m diameter. The required field-of-view is 
of about 10 X 10 arcmin^, with the goal of covering a four times 
larger area, so as to cover one entire rich cluster down to a relatively 
low redshift. The best band for tSZ observations will be centered 
on 150 GHz. At this frequency CCAT is expected to have Gaus- 
sian beam of 0.44 arcmin FWHM. The first step of noise setup is to 
convolve the maps with this beam. Then we add a Gaussian noise 
of 3^K/beam. This level of noise should be reached with about 6 
hours of exposure with CCAT. In the present study, we neglect in 
the tSZ maps any contamination, in particular we do not consider 
the presence of unresolved radio point sources. A detailed analy- 
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Figure 3. tSZ (upper panel) and X-ray (lower panel) maps of the ideal 
cluster. The side of the map corresponds to 2R^ir = 4 Mpc. In the bottom 
left comer of the tSZ map is the beam. Note that the scale is linear in the 
tSZ map and logarithmic in the X-ray one. 



sis of the contarn inations in the tSZ signal has been provided by 
iKnox et alj j2004) and bv lAghanim et al.l j2004h . 

As for X-ray observations, the Chandra satellite is currently 
providing imaging of superb quality, with a sub-arcsec resolution 
on axis. A proper simulation of X-ray observations should require 
generating spectra for each pixel, to be convolved with the response 
function of a given instrument. However, in order to apply our re- 
construction method we only need to generate X-ray surface bright- 
ness maps with a given number of events (photons), regardless of 
their energy. For this purpose we simulate the X~ray photon counts 
by using a Monte Carlo sampling of the surface brightness map. We 
fix to = 10* the total number of photons within the virial radius 
of the cluster, which is quite typical for medium-deep observations 
of relatively nearby clusters. Each photon event is generated in a 
particular pixel i, with probability 



Pi 



(10) 



where Si is the surface brightness of the pixel and the sum is ex- 
tended to all the pixels of the map. The number of expected counts 
into each pixel will be given by a Poisson probability distribution 
with mean rii = NPi . The conversion between counts and surface 
brightness is then given by E = ^ - Si /N, so that the total flux in 
the map is conserved. 

Clearly, this method of introducing noise in the X-ray maps 
only takes into account the statistical errors associated to finite 
exposures. However, it neglects the effects of any systematics 
(e.g., contribution of the instrumental and cosmic background, 
etc.) which should be included in a more realistic observational 
setup. A comprehensive description of the instrumental effects on 
the recovery of X-ray observables, calibrated on hydro dynamical 
simulations, has been provided by iGardini et all ( |2004) (see also 
iRasia et alj|200 ^. Probably the most serious limitation in our ap- 
proach is that we assume the absence of any background or, equiv- 
alently, that the background can be characterized and removed with 
arbitrary accuracy. 

In Figure[3]we show are the tSZ and X-ray images of the ide- 
alized cluster, once noise is added as described above. In Figure |2] 



we show the results of the deprojection of the noisy maps of the ide- 
alized cluster, for both the density and the temperature profiles. In 
order to estimate the errors in the deprojected profiles, we perform 
a MonteCarlo resampling of the projected X-ray and tSZ profiles: 
the value of the profile within each radial ring is randomly scattered 
according to a Gaussian distribution, whose width is given by the 
error associated to the noise introduced in the map. The la errors 
in the deprojected profiles is then obtained as the scatter within a 
set of 1000 deprojections of the MonteCarlo-resampled tSZ and 
X-ray profiles. 

The density is the best determined quantity, with uncertainty 
lower than 10 per cent. This is quite expected, owing to the sensitive 
dependence on gas density of both the X-ray signal (oc n^) and of 
the tSZ one (oc rie). The temperature has instead higher errors, of 
about 20-30 per cent. This is due to the fact that the X-ray signal 
has a weaker dependence on the temperature (only contained in 
the cooling function A(r)). For this reason, the determination of 
the temperature profile, independent of any X-ray spectroscopic 
analysis, is strictly related to the possibility of having high-quality 
tSZ data. 

The introduction of noise generates fluctuations in the tSZ and 
X-ray profiles which translate into variations of the recovered den- 
sity and temperature. Looking at the bottom panels of Fig. (2] pos- 
itive fluctuations in the density correspond to negative fluctuations 
in the temperature (and viceversa). Furthermore, any fluctuation in 
a given direction in one radial bin generally corresponds to a fluc- 
tuation in the opposite direction in an adjacent bin, within the same 
profile. This pattern in the fluctuations witnesses the presence of 
a significant covariance among nearby bins in the same profile and 
between the values of density and temperature recovered within the 
same radial bin. As for the covariance between neighbor bins, it is 
due to the onion-skin structure assumed in the deprojection: every 
time that a quantity is over(under)estimated in a radial bin, the de- 
projection forces the same quantity to be under(over)estimated in 
the adjacent inner bin, so as to generate the correct projected pro- 
file. As for the covariance between different profiles, it is mostly 
induced by the tSZ signal, which has the same dependence on both 
Tie and T. Although such oscillations are present for both density 
and temperature, they are smaller for the former, due to its faster 
decrease with radius. 



3.2 Maximum likelihood deprojection of the noisy maps 

We verified that using the maximum-likelihood technique, as de- 
scribed in Sect. 2.2, generally produces very similar results to those 
of the geometrical deprojection, at least when the regularization 
term, Creg is not included in the analysis. The results of this de- 
projection method on the polytropic /3-model are shown in Figure 
[4] where we also show the effect of introducing the regularization 
term. The effect of the regularization constraint is evident: most of 
the fluctuations, which are due to the degeneracies between fitting 
parameters, disappear and the deprojected profiles become much 
more regular, and with smaller errorbars in the profiles, while the 
accuracy of the reconstruction remains essentially unbiased. 

In order to study in detail the presence of correlations among 
the fitting parameters, we compute the correlation matrix, which is 
defined as dj — aij/aiaj, where aij is the covariance between 
the Xi and the Xj fitting parameters, while is the variance for the 
Xi parameter. The covariance matrix is computed along the Markov 
Chain. Therefore, dj is in our case a matrix with 30 x 30 entries. In 
Figures[5]and|6]we plot the entries of the correlation matrix for the 
density-density (DD), temperature-temperature (TT) and density- 
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Figure 4. Results of the maximum-likelihood deprojection on the density profile (left panel) and temperature profile (right panel) for the idealized cluster. In 
the upper part of each panel we show the correct profile (dashed curve) and the reconstructed profile with and without the regularization constraint (circles and 
triangles, respectively). In the lower part we show the fractional deviation of the recovered profiles from the true one. 
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Figure 5. The congelation matrix of density and temperature without the regularization constraint: density-density (left panel), temperature-temperature (central 
panel) and density-temperature (right panel). White pixels correspond to the presence of strong positive correlation, while black pixels are for strong anti- 
correlation. 
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Figure 6. The correlation matrix of density and temperature while using regularization constraint: density-density (left panel), temperature-temperature (central 
panel) and density temperature (right panel). The color-coding of the pixels is the same as in Figure|5] 
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temperature (DT) "blocks", before and after introducing the regu- 
larization term in tlie likelihood function, respectively. By defini- 
tion, the variance terms in the diagonal of the DD and TT matrices 
are characterized by the maximum correlation. On the contrary, the 
diagonal of the DT matrix has the maximum anticorrelation, thus 
demonstrating that any positive fluctuation in the recovered profile 
of one quantity corresponds to a negative fluctuation of the other 
quantity at the same radius. We also note that the next-to-diagonal 
terms in the DD and TT blocks have a degree on anticorrelation, 
thus explaining the fluctuating profile shown in Figs.|2]and|4] When 
the regularization is introduced, the correlations between density or 
temperature of adjacent bins is efficiently suppressed. 



4 APPLICATION TO SIMULATED CLUSTERS 

The sample of simulated galaxy clusters used in this paper has been 
extracted from the large-scale cosmological hydro-N-body simula- 
tion of a "concordance" ACDM model with Q,n = 0.3 for the 
matter density parameter at present time, Qa ~ 0.7 for the cos- 
mological constant term, Qh ~ 0.019 for the baryons den- 
sity parameter, h — 0.7 for the Hubble constant in units of 100 
km s'^Mpc^i and as = 0.8 for the r.m.s. density perturbation 
with in a top-hat sphere having comoving radius of Sh'^Mpc 
(see iBorgani et alj|2004l for further details). The run, performed 
with t he Tree+SPH code GADGET-2 jSpringel et al.ll200ll : ISpringell 
l2005h . follows the evolution of 480'^ dark matter particles and an 
equal number of gas particles in a periodic cube of size 192h~^ 
Mpc. The mass of the gas particles is mgas ~ 6.9 x W^h~^ AIq, 
and the Plummer-equivalent force softening is 7.5h~^ kpc at 
2 = 0. The simulation includes the treatment of radiative cool- 
ing, a uniform time-dependent UV background, a sub-resolution 
model for star formation and energy feedback from galactic winds 
jSpringel & Hemauisll2003h . The sample of clusters analyzed here 
includes 14 clusters extracted at jz = and having virial mass 
Mvir^ 4 X lO^^M©. For these clusters we compute the electron 
temperature, 

Te= ^ , (11) 

whose definition coincides with the mass-weighted one under the 
assumption of a fully ionized plasma. In the above equation, rui 
and ti are the mass and the temperature of the i-th gas particle. 

Due to the finite box size, the largest cluster found in the 
cosmological simulation has Te ~ 5 keV. In order to extend our 
analysis to more massive and hotter systems, which are mostly rel- 
evant for current tSZ observations, we include four more galaxy 
clusters having Mvir > lO^^^ fe-^Mp Ip and belonging to a differ- 
ent set of hydro-N-body simulations ('Borgan i et al.ll20()6h . Since 
these objects have been obtained by re-simulating, at high reso- 
lutio n, Lagrangian region s of a pre-existing cosmological simula- 
tion jYoshida et al]|200lh . they have a better mass resolution, with 
m-gas = 1.69 X 10^h~^ Mq, and a correspondingly smaller soften- 
ing of 5/i~^kpc at 2 = 0. These simulations have been performed 
with the same choice of the parameters defining star-formation and 

^ Here and in the following, the virial radius, Hvir, is defined as the ra- 
dius of a sphere centered on the local minimum of the potential, containing 
an average density, Pvir. equal to that predicted by the spherical collapse 
model. For the cosmology assumed in our simulations it is pvir — lOOpc, 
being pc the cosmic critical density. Accordingly, the virial mass, M^ir, is 
defined as the total mass contained within this sphere. 
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Table 1. Characteristics of the simulated clusters, for which the detail of the 
analysis are presented. Col. 1 : electron (mass-weighted) temperature; Col. 
3: virial mass; Col. 4: virial radius. 

feedback. The cosmological parameters also are the same, except 
for a higher power spectrum normalization, erg = 0.9. 

In the following, we will show detailed results for a subset of 
4 clusters. The basic characteristics of these four selected clusters 
are reported in Table [T] The first three of them are extracted from 
the cosmological hydrodynamical simulation, while C4 is one of 
the massive clusters simulated at higher resolution. C2, C3 and C4 
are typical examples of clusters at low, intermediate and high tem- 
perature, while CI is an interesting case to understand the effect of 
fore-background contaminations. We show in Figure |7] the X-ray 
surface brightness and Compton-j/ maps for these four clusters. All 
the maps are generated by placing the cluster at redshift z — 0.1, 
so that the maps, which extend out to r^ir, have an angular size 
ranging from about 9 arcmin for CI to 14 arcmin for C4. 

4.1 tSZ and X-ray maps 

Around each cluster we extract a spheric al region extending out 
to 6 Rvir- Following 'Diaferio et al. I COO?), we create maps of the 
relevant quantities along three orthogonal directions, extending out 
to about 2 Rvir from the cluster center, by using a regular 512x512 
grid. 

A number of different analyses, based on a joint deprojections 
of SZ and X-ray cluster maps under the assumption of axial sym- 
metry, indicate that the X-ray selection ten ds to favor objects which 
are elongated along the line-of-sight (e.g., iDe Filippis et alfcoOSl 
, and references therein). In order to control the effect of this selec- 
tion bias, we decided to choose the axes of projection to be aligned 
with the principal axes of inertia of the cluster. This will allow us to 
quantify the difference in the reconstructed profiles when the pro- 
jection direction is that of maximum cluster elongation. 

To derive these axes, we diagonalize the inertia tensor, which 
is given by 

N 

lij = ^{rirj)pl (12) 

p=0 

where i,j = 0,1,2 are the coordinate axes, is the i-th coordinate 
of the particle p with density pp and the sum is extended over all 
the gas particles. We weight each particle by pp so as to mimic the 
elongation in the X-ray emissivity. 

The eigenvectors of the I tensor provide the principal axes of 
the best-fitting ellipsoid. The semi-axes ai of this ellipsoid are pro- 
porti onal to square root of the corresponding eigenvalue ai oc y/\i 
(e.g. ,PlionisetalJll99lh . We choose the direction of projection z 
to be that corresponding to the largest semi-axis (i.e. the maximum 
elongation), while the y and x directions correspond to the medium 
and to the lower semi-axes, respectively. 

In the Tree-l-SPH code, each gas particle has a smoothing 
length hi and the thermodynamical quantities it carries are dis- 
tributed within the sphere of radius hi according to the compact 
kernel: 
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Figure 7. Maps of the X-ray surface brightness (left panel) and of the Comptonization parameter (right panel) for the CI to C4 simulated clusters (from top 
to bottom). Each map extends out to Ryir- Noise is added according to the recipe described in Section 3.1. Circles mark the regions which have been masked, 
due to the presence of detected substructures. 
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W{x) = 



-X)^ i < X ^ 1 

X>1 



(13) 



Here, itis x — r/hi with r the distance from the particle position. 
Therefore, we distribute the quantity of each particle on the grid 
points within the circle of radius hi centered on the particle. Specif- 
ically, we compute a generic quantity qjk on the (j, k) grid point 
as qjkdp = / q{r)dldp = ^ qi{mi / pi)Wi where dp is the pixel 
area, the sum runs over all the particles, and Wi oc / W{x)dl is the 



weight proportional to the fraction of the particle proper volume, 
rrii /pi, which contributes to the (j, k) grid point. For each particle, 
the weights are normalized to satisfy the relation ^ to^ = 1 
where the sum is over the grid points within the particle circle. 
When hi is so small that the circle contains no grid point, the par- 
ticle quantity is fully assigned to the closest grid point. 

As for the X-ray maps, they have been generated in the [0.5-2] 
keV energy band, by computi ng the emissivity of each gas particle 
with a Raymond-Smith code l lRavmond & Smithlll977h . assuming 
zero metallicity in the cooling function. 
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Figure 8. Application of tlie regularized maximum likelihood deprojection on the cluster CI, out to R^ir- The tree solid lines connecting dots with errorbars 
represent the reconstructed profile, for three orthogonal projection directions: along the x (squares), the y (triangles) and the z (circles) axes. En'orbars 
con'esponds the asymmetric 68 per cent confidence levels, computed from the distribution of values taken by the likelihood function along the Markov Chain. 
The dashed line represents the true 3-dimensional profile. The dotted line in the right panel shows the profile of the spectroscopic-like temperature. In the 
bottom panels, we plot the fractional deviation of the reconstructed profiles from the true electron temperature. 




Figure 9. The same as in Figure[8] but for the C2 cluster. 



Noise is finally added as described in Section [STI We fix the 
total number of photons in the virial radius to 10^ also for simulated 
clusters. For the tSZ map, we adopt a noise level of 10 /iK/beam 
for the objects having spectroscopic temperature Tsi > 4 K and 3 
/^K/beam for those having 3 K < Tsi < 4 K. 

4.2 Results 

Having tested the reliability of the deprojection method, with the 
regularization of the likelihood, we apply now this technique to 



the more realistic case of hydrodynamical simulations. In this case, 
a number of effects, such as deviations from spherical symmetry, 
presence of substructures and presence of fore/background con- 
taminating structures, are expected to degrade the capability of the 
deprojection to recover the three-dimensional profiles. 

We show in detail the results on the density and temperature 
deprojection the selected subset of four clusters presented in Table 
[T] while the whole set of 14 clusters will be used to assess on a 
statistical basis the efficiency with which the total gas mass can be 
recovered. 
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Figure 10. The same as in Figure[8] but for the C3 cluster. 
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Figure 11. The same as in Figure[8] but for the C4 cluster. 



The C2 and C4 objects are rather typical examples of our set 
of clusters. They are fairly relaxed and with a modest amount of 
substructures. As for the presence of substructures, they are well 
known to represent an important source of bias in the deprojec- 
tion, especially of the X-ray signal, which is highly sensitive to 
gas clumping. In order to remove this contaminating signal, we fol- 
low the same method that is often adopted in the analysis of ob- 
servational data. We first identify the detectable clumps by visual 
inspection of the X-ray maps. The corresponding regions are then 
masked out both in the X-ray and in the tSZ maps. The masked 
regions are excluded from the computation of the signals to be de- 
projected. This leads to an increase of the statistical uncertainties 
in those rings which have a significant overlap with the masked re- 



gions. Clearly, due to the finite photon statistics in the X-ray maps, 
small clumps may fall below the detection threshold, while their 
presence may still affect the emissivity. 

The recovered density and temperature profiles of C2 and C4 
are shown in Figures |9] and [TT] Once all the detectable clumps are 
masked out the reconstruction of the density profile is generally 
good, but with a systematic overestimate of ~5 per cent, that we at- 
tribute to a residual small-scale gas clumping. Although this effect 
is rather small, its presence highlights the need to have a sufficient 
photon-count statistics to identify gas inhomogeneities and remove 
their contribution in the deprojection procedure. The slight density 
overestimate corresponds, as expected, to a small underestimate of 
the temperature, which is forced by the requirement of reproducing 
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Figure 12. Projected surface gas mass density of the CI cluster along the 
z (top left), y (top right) and x (bottom left) directions, after removing the 
contribution from all gas particles which inside the virial region of the clus- 
ter. 



the tSZ signal, y oc TieTe- For these two objects we also note that 
there are rather small differences in the 3D profiles recovered from 
three orthogonal projection directions, thus indicating that they are 
almost spherical and without significant substructures along the dif- 
ferent projection directions. Errorbars are always of the order of a 
few percent, in both density and temperature. We stress that these 
very small errorbars, especially in temperature, are partly a conse- 
quence of the regularization constraint. 

As for the C3 cluster, we note that it has larger substructures 
which will have a stronger impact on the recovered profiles. Even 
after masking all the detectable substructures, we still have a num- 
ber of unresolved clumps. As expected, in this case the density 
profile (see Figure fTOt is overestimated by a larger factor, ~ 10 
cent, with a corresponding more significant underestimate of the 
temperature. The deviations of the deprojected profiles in the outer 
parts are also larger. This is due to stronger contaminations from the 
fore/background structures, which are both placed at the outskirts 
of the cluster and along the projection direction, in the cosmic web 
surrounding the cluster. In fact, both tSZ and X-ray maps are pro- 
duced by projecting a region of 6 Rvir in front and in the back of 
the cluster center. 

In this respect, the CI cluster is particularly interesting. Along 
the X-axis projection there is a merging groups along the same line 
of sight, at a distance of ~ 1.2Iiuir from the center of the main 
cluster. In Figure [T2] we show the projected mass surface density 
of the gas along the three projection directions, after removing the 
mass of the main cluster within R^ir- While the residual mass sur- 
face density is quite small along the x and y directions, a presence 
of a gas clump are shows up in the z projection. While this struc- 
ture provides a rather small contribution to the X-ray signal, its gas 
pressure is comparable to that of the main cluster, thereby signifi- 
cantly contaminating the tSZ effect signal. As a consequence, the 
density profile (see Figure [8j is essentially unaffected, while the 
temperature is clearly boosted by ~ 20 per cent with respect to 
the that obtained from the other two projections. Although this is 



a quite peculiar case, in which the secondary structure is relatively 
large and aligned with the main cluster along the line of sight, it 
illustrates the role of projection contamination from unidentified 
structures in recovering the 3D thermal structure of the ICM. 

We note in FigslMTTIthat the density profiles recovered from 
the projection of maximum elongation are overestimated at small 
radii, while they are underestimated in the outskirts. In order to 
quantify this effect, we show in the left panel of Figure[T3]the ratio 
between the true and the reconstructed density profiles, after aver- 
aging over the sample of simulated clusters. By averaging over all 
the projection directions, the density is generally overestimated by 
about 5 per cent at all cluster radii. This result is confirmed also by 
analyses performed on synthetic X-ray obse rvations of simulated 
clusters l lRasia et al.ll200a: iNagai et alj|2007h . On the other hand, 
the density profile reconstructed from the projection along the z 
axis is confirmed to be significantly larger than along the other 
direction in the very inner part, with an inversion at 0.2rvir- 
Indeed, the elongation causes the objects to appear more compact 
in the X-ray maps, which drive the density reconstruction. This 
boosts the deprojected central density, while depletes it in the out- 
skirts. 

As shown in the right panel of Fig. [13] the temperature is gen- 
erally underestimated by ^ 10 per cent out to ~ 0.7r„ir- At larger 
radii this underestimate increases, reaching a mean value of about 
20 per cent at r^ir, as a consequence of the relatively larger con- 
tamination by fore/background structures. The scatter is generally 
larger than the uncertainty introduced by the noise, so that it has 
to be considered as intrinsic to the measure. This scatter has dif- 
ferent origins, such as unresolved gas clumps, asphericity of the 
clusters, fore-background contaminations. In general, the tempera- 
ture recovered from the projection along the z axis is slightly larger 
than the one from the other two axes. The difference is more ap- 
parent in the central regions and becomes smaller in the outskirts. 
The reason for this behaviour is that the temperature reconstruction 
is more affected by the tSZ signal. Along the direction of maxi- 
mum elongation, this signal in enhanced since the ICM pressure 
is integrated along a larger path. The tSZ signal receives then an 
important contribution from cluster regions where the density is 
underestimated. As a result the reconstructed temperature is corre- 
spondingly increased to compensate for this effect. 

In Figs. [MTTI we also show the three-dimensional profile of 
the spect roscopic-like tempera ture (dotted line), which has been 
shown bv lMazzotta et al.l ( |2004|) to represent a quite close proxy to 
the actual X-ray temperature obtained from a spectroscopic fit (see 
alsO iVikhlinin 2006) . This temperature is computed as 



-1/2 



ET-iO — 3/2 



(14) 



with a ~ 0.75 and the sum extends over all the gas particles hav- 
ing internal energy larger than 0.5 keV. According to its definition, 
this temperature gives more weight to the low-temperature phase 
in a thermally complex ICM. This is the reason for the drop of the 
Tsi profile at the cluster center and for the wiggles which mark the 
positions of merging sub-clumps which are relatively colder than 
the ambient ICM. In general, the profile of Tsi are lower than those 
of the electron temperature, by an amount which is larger for hot- 
ter systems (see also Rasia et al. 2005). These figures highlight that 
the temperature profiles, as obtained from our deprojection analy- 
sis, are much closer to the mass-weighted temperature, which mea- 
sures the total thermal content of the ICM, than to T^i. An impor- 
tant consequence of this difference will clearly be the estimate of 
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Figure 13. The ratio between the reconstructed and the true (mass-weighted) temperature and density profiles. The shaded areas encompass the ztlcr regions 
of the recovered profiles over the ensemble of simulated clusters. The horizontally shaded area is for the projections along the z axis, while the vertically 
shaded area is for the projections along the other two axes. The black line shows the mean over all projections of all clusters. 
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Figure 14. The ratio between the recovered and the true values of the total 
gas mass for simulated clusters out to For each cluster we show the result of 
the deprojection along the three orthogonal directions, with the projection 
corresponding to the maximum elongation being marked with a filled circle. 
Errorbars con'espond to the Icr confidence level, by accounting for the full 
error correlation matrix when integrating the 3D gas density profiles. The 
horizontal dotted line shows the average value of the ratio. 



the total cluster mass from the application of hydrostatic equilib- 
rium. We will discuss the application of our deprojection method 
to cluster mass estimates in a forthcoming paper. 



4.3 Recovering the gas mass 

As a first application of our deprojection procedure, we compute 
the gas mass of the clusters, which is calculated simply by summing 
up the mass contained into each radial bin. Since the bins are not 
independent, the errors on the total gas mass have been calculated 
by using both variances and covariances of the values of the density 
at different radii. 

^ N N 

= iV^^IZl'^'"''^' ^^^^ 

1=0 j=0 

where am,ij is the covariance between the mass content of the i-th 
and of the j-th shells, directly obtained from the covariance be- 
tween the gas density in different bins (e.g., see left panel of Fig. 

Note that since the covariance between the density in adjacent 
bins is generally negative, neglecting it would lead to a systematic 
overestimate of the error on the mass. 

We give the results on the estimate of gas mass for the whole 
set of 14 simulated clusters. The small overestimate found for the 
density profiles is obviously propagated to the estimate of the total 
gas mass. The resulting bias turns out to be very small, and amounts 
to about 4 percent, with no obvious trend with the cluster mass. 
This demonstrate that residual gas clumping, after the removal of 
the substructures identified in the X-ray maps, has a small effect 
on the our capability of the recovering the total mass of the ICM. 
We note that cluster-by-cluster variance is often comparable to the 
"projection variance", i.e. to the differences found when projecting 
the same clusters along different directions. We also note that the 
uncertainties in the individual Mgas estimates, typically of the or- 
der of a few per cent, are smaller than the scatter. This indicates that 
the intrinsic scatter in the recovered gas mass is in fact associated 
to the deviations of the simulated clusters from perfect spherical 
symmetry. 
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Figure 15. The same of figure [T4l but as a function of cluster eUipticity (left panel), prolateness (central panel) and elongation along the line of sight (right 
panel). 



4.4 The effect of morphology 

In order to better understand how morphology affects the deprojec- 
tion, we plot in Figure [15] the recovered gas mass as a function of 
cluster ellipticity, prolateness and elongation along the line of sight. 
The ellipticity of a triaxial object is defined as: 



1 

2^ 



and the prolateness as: 
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The elongation is defined as the ratio between the semi-axis 
aligned with the line of sight and the larger of the other two semi- 
axes. 

The gas mass recovered from the projection along the princi- 
pal axis is generally lower than those from the other two projec- 
tions. This underestimate generally is anticorrelated with the elon- 
gation of the cluster, although with a substantial scatter From the 
left and the central panels of Figure [TS] we do not find any signif- 
icant correlation between the global 3D morphology of the clus- 
ters (prolateness and ellipticity) and the bias in the deprojection. 
Instead, as shown in the right panel, any effect in the gas mass re- 
covery is driven by the orientation of the cluster 



5 DISCUSSION AND CONCLUSIONS 

We have presented results of deprojection methods, aimed at re- 
covering the three-dimensional density and temperature profiles of 
galaxy clusters, by combining X-ray surface brightness and ther- 
mal SZ (tSZ) maps. The main aim of our analysis is to verify to 
what accuracy one can recover the thermal structure of the ICM by 
taking advantage of the different dependence of the X-ray and tSZ 
signal on the gas density and temperature, thereby avoiding per- 
forming X-ray spectroscopy. The two deprojection methods con- 
sidered are both based on assuming spherical symmetry of the clus- 
ters. 

The first one follows a geometrical approach, in which the 
3D profiles are recovered with an iterative procedure that depro- 
jects the observed images starting from the outermost ring and pro- 



ceeding inwards. The second method assumes the values of the 3D 
gas density and temperature profiles at different radii and computes 
from them the expected SZ and X-ray surface brightness which is 
then compared to the observations with a maximum likelihood ap- 
proach. In the computation of the likelihood, we also introduced 
a regularization term, which allows us to suppress spurious os- 
cillations in the recovered profiles. Using a Monte Carlo Markov 
Chain (MCMC) approach to optimize the sampling of the param- 
eter space, this second method also allows us to recover the full 
correlation matrix of the errors in the parameter fitting. 

The main results of our analysis can be summarized as fol- 
lows. 

• The application of both methods to an idealized spherical 
polytropic /3-model shows that the 3D profiles are always recov- 
ered with excellent precision (of about 3^ per cent), thus demon- 
strating that such methods do not suffer from any intrinsic bias. 

• The application of the maximum-likelihood method to hydro- 
dynamical simulations of galaxy clusters always provides depro- 
jected profiles of gas density and temperature, which are in good 
agreement with the true ones, out to the virial radius. We find a 
small (^ 10 per cent) systematic overestimate of the gas density, 
which is due to the presence of some residual gas clumping, which 
is not removed by masking out the obvious substructures identified 
in the X-ray maps. 

• The total gas mass is recovered with a small bias of 4 per cent, 
with a sizable scatter of about 5 per cent. This result shows that 
residual gas clumping should have a minor impact in the estimate 
of the total gas mass. We do not find any trend in the recovery of 
the gas mass with the total cluster mass. 

• The gas mass reconstructed along the maximum elongation 
axis is generally lower (by up to 10 per cent) with respect to the 
mass reconstructed along the other two projection axes, the size of 
this effect being larger for more elongated clusters. 

• The temperature is generally well recovered, with ~ 10 per 
cent deviations from the true one out to ~ 0.7-Rvir- The rather small 
size of this bias confirms that the combination with tSZ data is a 
valid alternative to X-ray spectroscopy for temperature measure- 
ments. The temperature reconstructed from the projection along the 
axis with maximum elongation is slightly higher than those from 
the other two axes, particularly in the inner regions. 

Our results confirm the great potentials of combining spatially 
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resolved tSZ and X-ray observations to recover the thermal struc- 
ture of the ICM. This approach has several advantages with respect 
to the traditional one based on X-ray spectroscopy. First, the tem- 
perature recovered from the fit of the X-ray spectra is known to 
provide a biased estimate of the total thermal content of the ICM, 
the size of this bias in creasing with the complexity of the p lasma 
thermal structure (e.g., iMazzotta et al.l2004l : rVikhlininl26o6h . Sec- 
ondly, X-ray surface brightness profiles can be obtained with good 
precision with a relatively small number of photon counts. Also, 
once the cosmic and instrumental backgrounds are under control, 
the surface brightness can be recovered over a large portion of the 
cluster virial regio ns, as already demonstrated with R OSAT-PSPC 
imagmg data (e.g., IVikhlininetalJl999l ; lNeumanrj20 05). Since the 
tSZ has the potential of covering a large range in gas density, then 
its combination with low-background X-ray imaging data will al- 
low one to recover the temperature profiles out to the cluster's out- 
skirts. 

A limitation of the analysis presented in this paper is that we 
did not include realistic backgrounds in the generation of the X- 
ray and tSZ maps. As we have just mentioned, there are reason- 
able perspectives for a good characterization of the X-ray back- 
ground. However, the situation may be more complicated for the 
tSZ background. In this case, conta minating signals from u nre- 
solved point-like radio sources (e.g., iBartlett & MeiiT]|2006l) and 
fore/background galaxy groups (e.g., iHallman et al. 2007) could 



affect the tSZ signal in the cluster outskirts. In this respect, the 
possibility of performing multi-frequency observations with good 
angular resolution will surely help in characterizing and removing 
these contaminations. 

Single-dish sub-millimetric telescopes of the next generation 
promises to provide tSZ images of clusters with a spatial resolution 
of few tens of arcsec, while covering fairly large field of views, 
with 10-20 arcmin aside, with a good sensitivity. At the same time, 
planned satellites for X-ray survevJTI will have the capability of 
surveying large areas of the sky with a good quality imaging and 
control of the background. These observational facilities will open 
the possibility of carrying out in survey mode high-quality tSZ and 
X-ray imaging for a large number of clusters. The application of 
deprojection methods, like those presented in this paper will pro- 
vide reliable determinations of the temperature profiles and, there- 
fore, to exploit their potentiality as tools for precision cosmology. 



ACKNOWLEDGMENTS. 

We would like to thank Stefano Ettori for his help with the geo- 
metrical deprojection algorithm. We also thank Antonaldo Diafe- 
rio, Sunil Golwala, Silvano Molendi, Manolis Plionis, Elena Rasia 
and Paolo Tozzi for useful discussions. This work has been par- 
tially supported by the PD-51 INFN grant and by the ASI-INAF 
1/023/05/0 grant. E.P is an ADVANCE fellow (NSF grant AST- 
0649899). 



REFERENCES 

Aghanim N., Hansen S. H., Lagache G., 2004, ArXiv Astro- 
physics e-prints 



eROSITA: http : / / www . mpe . mpg .de /projects . html#erosita 
* EDGE: |http ://ibis.rm.iasf.cnr. it /EdgeOverview . htm| 



Ameglio S., Borgani S., Diaferio A., Dolag K., 2006, MNRAS, 
369, 1459 

Bartlett J. G., Melin J.-B., 2006, A&A, 447, 405 
Birkinshaw M., 1999, Phys. Rep., 310, 97 

Bonamente M., Joy M. K., LaRoque S. J., Carlstrom J. E., Reese 
E. D., Dawson K. S., 2006, ApJ, 647, 25 

Borgani S., 2006, ArXiv Astrophysics e-prints 

Borgani S., Dolag K., Murante G., Cheng L.-M., Springel V., 
Diaferio A., Moscardini L., Tormen G., Tornatore L., Tozzi P., 
2006, MNRAS, pp 270-+ 

Borgani S., Murante G., Springel V., Diaferio A., Dolag K., 
Moscardini L., Tormen G., Tornatore L., Tozzi P., 2004, MN- 
RAS, 348, 1078 

BouchetL., 1995, A& AS, 113, 167 

Buote D. A., 2000, ApJ, 539, 172 

Carlstrom J. E., Holder G. R, Reese E. D., 2002, ARAA, 40, 643 
Cavaliere A., Fusco-Femiano R., 1976, A&A, 49, 137 
Cavaliere A., Lapi A., 2006, ApJ, 647, L5 

Croston J. H., Arnaud M., Pointecouteau E., Pratt G. W., 2006, 

A&A, 459, 1007 
De Filippis E., Sereno M., Bautz M. W., Longo G., 2005, ApJ, 

625, 108 

Diaferio A., Borgani S., Moscardini L., Murante G., Dolag K., 
Springel V., Tormen G., Tornatore L., Tozzi P, 2005, MNRAS, 
356, 1477 

Dore O., Bouchet F. R., Mellier Y., Teyssier R., 2001, A&A, 375, 
14 

Ettori S., De Grandi S., Molendi S., 2002, A&A, 391, 841 
Gardini A., Rasia E., Mazzotta P., Tormen G., De Grandi S., 

Moscardini L., 2004, MNRAS, 351, 505 
Gelman A., Rubin D., 1992, Stat. Science, 7, 457 
Gilks W., Richardson S., Spiegelhalter D., 1996, Markov Chain 

Monte Carlo in practice. Chapman and Hall 
Hallman E. J., O'Shea B. W., Burns J. O., Norman M. L., Hark- 

ness R., Wagner R., 2007, ArXiv e-prints, 704 
Hastings W., 1970, Biometrika, 57, 97 
Knox L., Holder G. P, Church S. E., 2004, ApJ, 612, 96 
Kotov O., Vikhhnin A., 2006, ApJ, 641, 752 
Kriss G. A., Cioffi D. F, Canizares C. R., 1983, ApJ, 272, 439 
LaRoque S. J., Bonamente M., Carlstrom J. E., Joy M. K., Nagai 

D., Reese E. D., Dawson K. S., 2006, ApJ, 652, 917 
Lee J., Suto Y, 2004, ApJ, 601, 599 
Lewis A., Bridle S., 2002, Phys. Rev. D, 66, 103511 
MacKay D., 1996, Markov Chain Monte Carlo in practice. Chap- 
man and Hall 

Mazzotta P., Rasia E., Moscardini L., Tormen G., 2004, MNRAS, 
354, 10 

McLaugWin D. E., 1999, AJ, 117, 2398 

Metropolis N., Rosenbluth A., Rosenbluth M., Teller A., Teller E., 

1953, Journal of Chemical Physics, 21, 1087 
Morandi A., Ettori S., Moscardini L., 2007, MNRAS, pp 541-1- 
Nagai D., Vikhlinin A., Kravtsov A. V., 2007, ApJ, 655, 98 
Neal R. M., , 1993, Probabilistic inference using Markov Chain 

Monte Carlo methods 
Neumann D. M., 2005, A&A, 439, 465 

Piffaretti R., Jetzer P, Kaastra J. S., Tamura T., 2005, A&A, 433, 
101 

Plionis M., Barrow J. D., Frenk C. S., 1991, MNRAS, 249, 662 
Pointecouteau E., Hattori M., Neumann D., Komatsu E., Matsuo 

H., Kuno N., Bohringer H., 2002, A&A, 387, 56 
Pratt G. W., Arnaud M., 2005, A&A, 429, 791 
Puchwein E., Bartelmann M., 2006, A&A, 455, 791 



Raftery 2003, Information Theory, Inference, and 
Learning Algorithms. Cambridge University Press 
(http://www.inference.phy.cam.ac.uk/mackay/itprnn/book.htmI 

Rasia E., Ettori S., Moscardini L., Mazzotta R, Borgani S., Dolag 
K., Tormen G., Cheng L. M., Diaferio A., 2006, MNRAS, 369, 
2013 

Rasia E., Mazzotta P., Borgani S., Moscardini L., Dolag K., Tor- 
men G., Diaferio A., Murante G., 2005, ApJ, 618, LI 
Raymond J. C, Smith B. W., 1977, ApJS, 35, 419 
Reblinsky K., 2000, A&A, 364, 377 
Rosati R, Borgani S., Norman C, 2002, ARAA, 40, 539 
Sebring T. A., Giovanelli R., Radford S., Zmuidzinas J., 2006, in 
Ground-based and Airborne Telescopes. Edited by Stepp, Larry 
M.. Proceedings of the SPIE, Volume 6267, pp. 62672C (2006). 
Cornell Caltech Atacama Telescope (CCAT): a 25-m aperture 
telescope above 5000-m altitude 
Sereno M., 2007, ArXiv e-prints, 707 
Springel V., 2005, ArXiv Astrophysics e-prints 
Springel V., Hernquist L., 2003, MNRAS, 339, 289 
Springel V., Yoshida N., White S., 2001, New Astronomy, 6, 79 
Sunyaev R. A., Zeldovich Y. B., 1972, Corrmients on Astrophysics 

and Space Physics, 4, 173 
Vikhhnin A., 2006, ApJ, 640, 710 
Vikhlinin A., Forman W., Jones C, 1999, ApJ, 525, 47 
Vikhlinin A., Markevitch M., Murray S. S., Jones C, Forman W., 

Van Speybroeck L., 2005, ApJ, 628, 655 
Voit G. M., 2005, Reviews of Modem Physics, 77, 207 
Yoshida N., Sheth R., Diaferio A., 2001, MNRAS, 328, 669 
Zaroubi S., Squires G., de Gasperis G., Evrard A. E., Hoffman Y, 

SilkJ., 2001, ApJ, 561,600 
Zaroubi S., Squires G., Hoffman Y, Silk J., 1998, ApJ, 500, L87+ 
Zhang T-J., Wu X.-R, 2000, ApJ, 545, 141 



Joint deprojection of SZ and X-ray images 



- ^.o.Mr [Me] 

5x1014 1015 



0.4 


1 III 


1 1 


1 

O 


1 


1 


0.3 


o o 


o 


o 






0.2 


o 

o 


o 


o 

o 


o 


— 


0.1 


1 III 


1 1 1 


1 


1 


1 




5x1013 


1014 

"^gas.vir 


[Mo] 







